Time–space fractional governing equations of transient groundwater flow in confined aquifers: Numerical investigation          下载免费PDF全文

Tongbi Tu  Ali Ercan  M. Levent Kavvas 《水文研究》2018,32(10):1406-1419

In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field.  相似文献   

黑子磁流体静力学的一个模型     

刘庆忠  宋慕陶 《天文学报》1996,37(1):35-42

本文对Ｏｓｈｅｒｏｖｉｃｈ的黑子返回磁通量模型作了适当的修改，使用黑子中心作为边界条件，用五种观测结果，导出了理论模型所需的五个主要参量，用半经验方法求得了黑子静力学模型的磁场、压力和温度等物理量．将此模型应用于一个中等大小的圆形对称黑子，可得到一个特解，结果发现我们的模型既能满足黑子的磁性质，同时又能满足合理的热力学量分布．  相似文献   

有关边界层问题中一类奇异积分方程新的存在性结果     

石丽莉  张高勋 《成都信息工程学院学报》2006,21(6):891-896

对边界层理论新结果中出现的一类奇异积分方程w（t）=∫1 t（1-s）（λ＋λs＋s）/w（s）ds＋（1-t）∫t 0s/w（s）ds,t∈（0,1）进行讨论，并得出了上述方程在λ∈（-1/2,0）上正解存在性的新结果。  相似文献   

关于Diophantine方程的解     

乐茂华 《广东海洋大学学报》2004,24(4):47-48

对于正整数a ,设S(a)是Smarandache函数。证明了 :方程S(1·2 ) +S(2·3) +… +S(x(x +1) ) =S(x(x +1) (x +2 ) /3)仅有正整数解x =1。  相似文献   

Modelling of deformation response and chemo‐mechanical coupling in chalk     

S. Pietruszczak  D. Lydzba  J. F. Shao 《国际地质力学数值与分析法杂志》2006,30(10):997-1018

A constitutive relation is derived for describing the mechanical response of chalk. The approach is based on a phenomenological framework which employs chemo‐plasticity. The properties of the material are assumed to be affected by the physico‐chemical processes that occur through the interaction between the skeleton and the pore fluid. The underlying mechanism is discussed by invoking a micromechanical analysis. The performance of the framework is illustrated by examining the evolution of mechanical characteristics in the presence of different pore fluids. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

CCD联合平差方法及其初步实例检验     

YU Yong LI Jinling TANG Zhenghong WANG Guangli ZHAO Ming 《中国科学院上海天文台年刊》2003,(1)

在常规照相天体测量工作中 ,对暗天体的定位是采用逐级定标星过渡的方法实现的。由于星表系统差和局部差的存在 ,最终定位结果可能包含难以把握的误差积累。在CCD小视场观测情况下时常难以找到足够数量的定标星。鉴于此 ,推导了CCD观测联合平差方法的严格矢量表达式 ,以期通过共同星的过渡 ,实现大天区统一平差。原理上此方法可实现相对于河外天体的定位 ,以削弱星表系统差和局部差对定位结果的影响。讨论了CCD联合平差方法实现中的问题 ,包括底片常数模型选择、切点位置改正、共同星的较差改正和法方程解算方法等 ,并分析了具体的处理措施。用实例初步检验了CCD联合平差方法的效果 ,表明通过对多张CCD观测的联合平差 ,可以削弱单张CCD观测参考星数目过少、单张CCD观测参考星分布不均、个别参考星存在位置偏差和局部区域参考星存在系统性位置偏差等不利因素对归算结果的影响 ,进而达到扩大视场和提高归算精度的目的  相似文献   

褶皱岩层初始挠曲的分叉特性     

杨成永  陈至达  杨锡禄 《大地测量与地球动力学》1995,(3)

基于拖带坐标法建立了岩层大变形的一维非线性微分方程，然后用解析与数值方法研究了变刚度、可伸缩岩层在侧向压力及自重作用下的分叉特性，初步揭示了褶皱岩层初始挠曲的形成机理。  相似文献   

凤太矿田层控铅锌(铜)矿床的成矿物质来源探讨     

隗合明  吕仁生 《地球科学与环境学报》1990,(1)

秦岭凤太矿田层控铅锌(铜)矿床的金属物质、硫和成矿溶液主要来自深部基底的岩石,属海底喷流—沉积成因矿床。从这一理论认识出发可进一步分析控矿地质条件及今后在该区有效地寻找同类矿床。  相似文献   

Restricted Three-Body Problem: An Approximation of its General Solution – Part One – the Manifold of Symmetric Periodic Solutions     

Kostantinos Papadakis  Constantine Goudas 《Astrophysics and Space Science》2006,305(2):99-124

This paper gives the results of a programme attempting to exploit ‘la seule bréche’ (Poincaré, 1892, p. 82) of non-integrable systems, namely to develop an approximate general solution for the three out of its four component-solutions of the planar restricted three-body problem. This is accomplished by computing a large number of families of ‘solutions précieuses’ (periodic solutions) covering densely the space of initial conditions of this problem. More specifically, we calculated numerically and only for μ = 0.4, all families of symmetric periodic solutions (1st component of the general solution) existing in the domain D:(x ₀ ∊ [−2,2],C ∊ [−2,5]) of the (x ₀, C) space and consisting of symmetric solutions re-entering after 1 up to 50 revolutions (see graph in Fig. 4). Then we tested the parts of the domain D that is void of such families and established that they belong to the category of escape motions (2nd component of the general solution). The approximation of the 3rd component (asymmetric solutions) we shall present in a future publication. The 4th component of the general solution of the problem, namely the one consisting of the bounded non-periodic solutions, is considered as approximated by those of the 1st or the 2nd component on account of the `Last Geometric Theorem of Poincaré' (Birkhoff, 1913). The results obtained provoked interest to repeat the same work inside the larger closed domain D:(x ₀ ∊ [−6,2], C ∊ [−5,5]) and the results are presented in Fig. 15. A test run of the programme developed led to reproduction of the results presented by Hénon (1965) with better accuracy and many additional families not included in the sited paper. Pointer directions construed from the main body of results led to the definition of useful concepts of the basic family of order n, n = 1, 2,… and the completeness criterion of the solution inside a compact sub-domain of the (x ₀, C) space. The same results inspired the ‘partition theorem’, which conjectures the possibility of partitioning an initial conditions domain D into a finite set of sub-domains D i that fulfill the completeness criterion and allow complete approximation of the general solution of this problem by computing a relatively small number of family curves. The numerical results of this project include a large number of families that were computed in detail covering their natural termination, the morphology, and stability of their member solutions. Zooming into sub-domains of D permitted clear presentation of the families of symmetric solutions contained in them. Such zooming was made for various values of the parameter N, which defines the re-entrance revolutions number, which was selected to be from 50 to 500. The areas generating escape solutions have being investigated. In Appendix A we present families of symmetric solutions terminating at asymptotic solutions, and in Appendix B the morphology of large period symmetric solutions though examples of orbits that re-enter after from 8 to 500 revolutions. The paper concludes that approximations of the general solution of the planar restricted problem is possible and presents such approximations, only for some sub-domains that fulfill the completeness criterion, on the basis of sufficiently large number of families.  相似文献   

排序对模糊度解算中降相关性能的影响分析     

鲁铁定  汪鑫  卢立果  徐兆祥 《大地测量与地球动力学》2020,40(5):470-475

针对不同排序算法对模糊度解算存在降相关性能影响的问题,从理论上分析了自然升序法、对称旋转法及扰动升序法的降相关原理,并基于模拟数据和实测数据,从降相关时间、搜索时间、总体耗时、Bootstrapping成功率及条件数5个方面对3种算法进行对比分析。结果表明,降相关效率与搜索椭球压缩程度呈负相关关系,搜索椭球压缩程度越高,降相关效率越低;对于不同的排序算法,提高降相关性能的关键在于减少降相关时间及对条件方差按一定方向排序,进而提高搜索效率。  相似文献